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Abstract. Phenomenological and theoretical aspects of fragmentation for elemen- 
tary particles (resp. nuclei) are discussed. It is shown that some concepts of classical 
fragmentation remain relevant in a microscopic framework, exhibiting non-trivial prop- 
er) | erties of quantum relativistic field theory (resp. lattice percolation). 
OV 
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' | ■ Introduction 

At first sight, fragmentation appears to be a typical classical process [1], by contrast 
with a quantum (or microscopic) one. You take a "big" object of size A, you break it 
into pieces of different sizes and - if you are a scientist - you count the mean number 
Nj of pieces of size j after a certain number of events. In rather simple but general 
cases [2] , one may introduce a time-dependent description in terms of successive binary 
splittings which can be formulated in terms of a "gain-loss" equation, namely: 
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k=j+l 

where Wjk is the binary fragmentation weight of an intermediate fragment of size k 
into j and k — j and, 

1 J_1 

£=1 

represents the total loss rate for the fragment species of size j. Indeed, there exist 
well-known papers [3] which discuss various mathematical solutions of equations of the 
type (1-2), with applications to e.g. depolymerisation, breaking of clusters etc... 



1 



For fragmenting quantum objects, such as elementary particles or nuclei, however, 
classical fragmentation concepts and equations are not necessarily relevant. In these 
cases, the fragmenting "big object" and the parameter A correspond to a particle jet 
and its total energy before its fragmentation into pions and other particles. For nuclei, 
it is an excited nucleus of atomic mass A fragmenting into smaller ones (including 
individual nucleons or a particles etc.). Note that in the former case, one has to 
replace the summation by an integral in equations (1-2), as also considered in Refs.[3]. 

There are some basic obstacles in front of us if we want to discuss quantum frag- 
mentation in terms of classical concepts. For instance, for elementary particles and 
nuclei, the excited quantum state which characterizes the system before its fragmen- 
tation is governed by the same interactions (called virtual) than those (called reai) 
responsible for its subsequent fragmentation. It is thus by no means obvious that the 
classical fragmentation structure could emerge from the quantum environment of the 
process. For instance, in field theory, the quantum fluctuations (e.g. loop Feynman 
diagrams) and more generally the renormalization procedure (necessary to give a real- 
istic meaning to the perturbative calculations) do not necessarily lead to an equation 
like (1). Quite unexpectedly, equation (1) is however useful for Quantum Chromody- 
namics (QCD), the theory of the fundamental strong interactions between quarks and 
gluons, as we shall see further on. 

Our aim in this contribution is to give a brief survey of how one can cast a bridge 
between quantum and classical concepts of fragmentation in this context. We will 
briefly analyze two specific cases where the phenomenological analysis can be sup- 
ported by a theoretical model. In section 1, Quark jet fragmentation will be analyzed 
in the framework of Quantum Chromodynamics (QCD). In section 2, Nuclear multi- 
fragmentation will be discussed in the framework of 3-dimensional lattice percolation. 



1. Quark Jet Fragmentation 



The best occasion where we can observe and measure quark jet fragmentation is the 
e + e~ annihilation into hadrons (mainly pions and their decay products) at high energy 
when an intermediate Z° boson is formed, thereafter decaying into quark-antiquark 
(qq) pairs. Such experiments have been performed at the LEP accelerator at CERN. 

During such reactions, in a first step lasting less than 10 -24 second, an interme- 
diary Zq boson is formed and decays into a qq pair, often followed by the subsequent 
formation of a third (gluon) jet. These jets form many gluons and other qq pairs. This 
stage is well described by QCD calculations with a small effective coupling constant [4] 
and can thus be studied in a quantitative theoretical framework. Then, in a later 
stage of the fragmentation process, quarks and gluons recombine into hadrons in an 
unknown way, only described by modelization. The deep fundamental reason of this 
is that as, the effective coupling constant of QCD, happens to be time-dependent as 
a consequence of quantum fluctuations. One has 
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where 6, Aqcd are fundamental constants and Q ~ 1/Time. One says that the coupling 
constant is "running" in this theory, beeing small at short times and becoming large 
later on. One also speaks of "asymptotic freedom" when time is short and "infrared 
slavery" at long times, since the elementary quanta, q, q, g, (called partons) are quasi 
free at short times and become tightly bound at long times and confined into hadrons. 
In this limit, the theory is in a non-perturbative regime and its complete solution is 
not yet known. 

At short times, the predominance of quantum fluctuations and virtual interaction 
effects makes difficult a classical fragmentation picture of a quark jet. However such a 
description emerges from the calculations after using a set of non-trivial properties of 
QCD, the field theory of Gauge Fields describing strong interactions between partons. 
This theory possesses an internal symmetry group, local in space-time, which is the 
"color" group SU3. We have no place for giving the full derivation of the gain- loss 
equation associated with QCD, but let us only describe it. The main issue lies in a 
system of equations which takes the form of a continuous "gain-loss" expression similar 
to (1,2). 

One writes: 
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where D A is the probability distribution of finding a "quantum" B or parton (g, q or 
q) in the fragmentation of the initial "quantum" A (g, q or q), with the fraction z of its 
total momentum. Note that Eq. (4) can be obtained, after some manipulation, from 
a continuum limit of (1) by choosing: 

-N^DU/A); W 3k = \p{j/k)- t = C 
J K 

Following the analogy with the classical process obeying Equations (1),(2), one may 
interpret the first term in the integral as a "gain" term where the parton B is obtained 
via first fragmentation of A into an intermediate parton C. Standard QCD calculations 
give a specific prediction for the weights P A (x) and thus for the solution of equation 
(4) . Note that the theory also leads to a precise re-definition of the "evolution" variable 
£, namely 

/" Qmax as(Q) dQ 1 ln(Q max /A QCD ) ^ 
* Jq 2tt Q 2wb n in (Q/Aqcd) ' 1 ' 

From various experimental analyses it has been shown that equation (4) gives 
a good description of the energy spectrum of jet fragmentation. However a model- 
dependent piece of information has to be added since one measures hadrons and not 
partons in the final state. In other physical configurations (structure functions instead 
of fragmentation) the same equation holds and can be tested with great success [4]. 

Note the interesting property of equation (4) that it can be exactly solved by 
the method of moments (or Mellin-transform) . Introducing the moments [M g ]^ = 
D A (z), one gets in matrix form: 

[M,]=exp , (6) 



fo dz 
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where the matrix elements [H q ] A are the q— moments of the weights P A . 

The emergence of a tree structure in jet fragmentation is not only based on the 
energy spectrum given by equation (4). Many other observables lead to the same struc- 
ture, while the detailed analyses show that it always implies a non-trivial property of 
both the quantum and group invariance properties of the theory. As an illustration, 
it was recently shown [5] that the multiplicity fluctuations associated with jet frag- 
mentation possess a dynamical structure, similar to the intermittency phenomena in 
hydrodynamics, which was predicted in particle physics some time ago [6]. This brings 
an interesting analogy between quantum fragmentation and intermittent fragmentation 
models, which appear in various classical or non classical systems such as spin-glasses, 
polymer diffusion, multi-particle production [7]. 

Note, however, that "differential" fragmentation observables beyond the energy 
spectrum and the functions D A (z) are more dependent on the unknown "hadroniza- 
tion" phase of partons and thus at present more model-dependent. The research is 
going on in this field. 

2. Nucleus Multifragmentation 

The physical understanding of nuclear multifragmetation is much less advanced than in 
QCD jet fragmentation. Experimentally , it is only recently that systematic data on the 
decay products of fragmented nuclei hence become available thanks to An— detectors 
at nuclear accelerators [8] . Even then, the difficulty remains of specifying without am- 
biguity the excited system which multifragments, and separating its fragments from 
the pre-equilibrium particles. Theoretically, the kinematical conditions of nuclear mul- 
tifragmentation, e.g. the incident energy and the multicomponents of the final state, 
are far from a known regime of nuclear forces. One has thus to rely on models, which 
are useful for the experimental investigation and may lead to a deeper understanding 
of the multifragmentation phenomenon. 

Among the proposed models, let us choose and discuss the one based on 3- 
dimensional percolation on a finite lattice[9,10]. Due to its particular simplicity, though 
unexplained on a purely nuclear-theoretical framework, it will allow us to develop on 
the links between a microscopic description and classical fragmentation concepts. This 
part of the talk comes from a recent study done in collaboration with Bertrand Giraud 
in Saclay[ll]. 

In the 3— d percolation model, the excited nuclear system of mass A (in nucleon 
mass units) is modelized by a finite lattice of volume A with nucleon on sites and 
nearest-neighbour bonds. Multifragmentation results from a breaking of these bonds 
in proportion of the energy release in the system by the reaction. Event-by-event some 
bonds are randomly preserved, corresponding in average to a ratio p, (0 < p < 1) of 
all bonds, while a certain number of fragments are formed, giving rise to a statistical 
distribution of fragments as a function of their size i, (1 < i < A). This distribution is 
in good agreement with recent data, if one replaces the unknown parameter p by an 
observable input for each event, e.g. the multiplicity of fragments. Another parameter 
is introduced corresponding to a site occupation probability, but we will stick to the 
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original model for simplicity. 

Our aim[ll,12] is to answer the following question: is it possible to find a classical 
time-dependent description of multifragmentation which would give, at least within 
some approximation, the same prediction than 3— d percolation concerning the distri- 
bution of fragments? In some sense, we are looking for an eventual restoration of the 
time variable in the percolation problem where such a reference scale is absent. More- 
over, the question behind this is whether percolation could be described by a linear set 
of equations similar to Eqns. (1-2). If such is the case, one could look for new scaling 
laws, in much the same way as in the case of jet fragmentation where they correspond 
to scaling properties of QCD. 

Technically, our workfll] starts with the quest of a general, albeit approximate, 
solution of the gain-loss equations (1-2). The idea is to exhibit properties which would 
be independent of the particular choice of weights Wjk, (which we do not yet know for 
percolation) . Then, one looks for the same properties in some range of the percolation 
model to test an eventual compatibility. Our conclusion is that indeed this compat- 
ibility can be achieved during the short time evolution of the system. Let us sketch 
how this can be proven. 

As any such linear system, the solution of Eqns (1-2), is otained via exponentiation 
once the eigenvalues and eigenvectors of the matrix [Wjk — Cjdjk] are determined. Note 
that the matrix is triangular, and thus the elements Cj on the diagonale are the exact 
eigenvalues. However the eigenvectors are unknown, except that they form also a 
triangular matrix. In fact, we were able to prove that these vectors take the quite 
general form of eigenmoments, namely they are of the form M q ^, where M q is the 
moment of rank q of the mean distribution of fragments and q(j) > l,j = 1, A are 
particular, but not necessarily integer, values of the rank. One has: 

q(l) = Kq(2)<q(3)...<q(A). 

This property is obtained by inspection of large matrices for which the problem is 
very similar to the QCD case (see Eq. (4-6)) where the moments M q give an exact 
diagonalization in the space of particle momentum. What was verified also for matrices 
of limited size, is that the diagonalization by moments remains true for a discretized 
set of values of the rank, a set q(j). Note however the model dependence of the set 
q(j) except for q(l) = 1 which is dictated by mass conservation. A limitation of the 
method was found in the case of the so-called "shattering transition", see [3], which 
needs a special treatment. With these limitations in mind, the eigenmoment property 
is general enough to be tested, e.g. in the case of percolation. 

For this purpose, we remark that, if they are identified as eigenmoments, the M q s 
are linked by linear relations in Log-Log plots, and their explicit time dependence 
disappears. We are thus led to display in the same way the moments obtained from 
the percolation model, choosing for instance Mi for reference, see the figure. Different 
moments are displayed (with q = 1, 1.5, 2, 3, 4, 5) and show the interesting feature of a 
quasi-linear dependence for the values q = 3, 4, 5, given the fact that for q = 1 (mass 
conservation) and q = 2 (reference scale) the linear relation is trivial. It is clear from 
this figure that a quasi-linear form is obtained between p = 1 and p = p c , where p c is 
the critical value above which, in the continuous limit, an infinite percolation cluster 
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is formed. Indeed, the figure shows the dominant contribution of the cluster of largest 
mass to the averaged moments. This largest cluster is, for finite size problems, the 
representative of the infinite cluster when p > p c . 

Notice that the moments implied by the rate equations are the full moments, 
including the largest fragment, while in usual analyses of percolation models[9,10], 
scaling properties are investigated with moments modified by the subtraction of the 
largest cluster. Moreover, in such traditional analyses of percolation, the reference 
time scale is generally given by the moment Mq or a similar variable related to the 
multiplicity of fragments. The comparison and compatibility of our approach with 
such analyses is an open problem of some interest. 

Conclusion and prospects 

The problem we want to discuss in conclusion of the study of the particle and 
nuclei fragmentation is whether classical concepts of fragmentation could serve as an 
unifying phenomenological picture at the microscopical level where quantum states 
and fields are involved. For this sake, let us discuss the striking common features and 
differences between the two examples we have treated. 

In the case of elementary particles and field theory, it is known that only the 
short time evolution of a jet is accessible to perturbative calculations. More precisely, 
it is only the time derivative of the fragmentation functions (or momentum distribu- 
tions of partons) which can be calculated exactly at first order in the quantum loop 
expansion. At any given fixed time, however, the knowledge of these functions de- 
pend on an expansion at all orders (for which the renormalization group properties 
can be invoked[4]). At long times it involves the unknown transformation of partons 
into hadrons. In other terms, the change in the vacuum structure (from partons to 
hadrons) prevents one from a complete theoretical understanding of the fragmentation 
process. 

In the case of nuclear multifragmentation, one does not possess a comparable the- 
oretical framework. However, the indications from the 3— d percolation model shows 
that there could be a similar property at short times, namely the possibility of a mul- 
tilinear evolution of fragmentation. This time range corresponds to the situation when 
the fragmenting force is mild enough to preserve the existence of at least one large 
cluster, that is when p < p c . What seems to be remarkable is the complementarity of 
percolation with the previous case. Percolation is well determined in a given region of 
the parameter, namely near the critical percolation value p c . In other terms, fragmen- 
tation is better determined at a given "time", while its "time-dependence", e.g. the 
relation between p— and time-evolution is not trivial. This is just the opposite of the 
field theoretical case. It is tempting to confront the methods used in the two cases in 
such a way that the stated complementarity could hopefully be used to explore the 
shadow regions of both processes. 

In conclusion, if fragmentation concepts could acquire some kind of universality, 
one can hope to find new methods to overcome the difficulties of the physical description 
of fragmentation in the quantum world. Some property of nuclear multifragmentation 
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could be useful for the hadronization problem of parton jets, as well as perturbative 
methods of Quantum Field Theory could help solving some nuclear fragmentation 
puzzles. More work in these directions is deserved. 
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Percolation analyzed with the M 2 time scale 



Relative strengths of moments M q , q = 0, 1, 1.5, 2, 3, 4, 5, as functions of M 2 in a Log- 
Log plot. Data taken from 3-d bond percolation on a 6*6*6 lattice. The corresponding 
values of the bond survival probability p are shown on the horizontal axis. Its critical 
value is p c = .25. Full lines: moments. Dashed lines: contributions of the largest 
cluster. Dashed-dotted line: the reference moment M 2 . Notice that a linear behaviour 
is approximately obtained for < p < p c and q = 3,4, 5 (Figure from [11]). 
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